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Abstract 

Using a gauge covariant operator technique we deduce the path integral for a 
charged particle in an arbitrary stationary magnetic held, verifying the "midpoint 
rule" for the discrete form of the interaction term with the vector potential. 
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The path integral for a charged particle in a magnetic field gives us a remarkable 
example of the subtleties of the discretized time approach and the relation between 
path integrals and Brownian motion. When one starts with the continuum definition 
of the path integral and tries to perform the functional integration by time slicing one 
is faced with the discretization of the term / dt x • A, that is the interaction with the 
vector potential A(x), and that for each infinitesimal segment x n+1 — x n is given by 
(x n+ i — x n )-A(x). Naively one could think that it would not matter whether A(x) 
is evaluated at x n , x n+ i or somewhere in between. But it happens that for the path 
integral to satisfy the Schrodinger equation we must take the average of A(x) between 
x ra and x n+ i ("midpoint rule"). If we disregard this rule we are forced to introduce an 
extra term in the Schrodinger equation, breaking the gauge invariance. The reason for 
that is in the Brownian character of the paths in the functional integral Jl]. Another 
approach to deduce the discretized path integral is to start from the Hamiltonian 
formulation of quantum mechanics and define it as an infinite product of operators 
where careful use of the Trotter formula for a proper "x — p" ordering leads to the 
"midpoint rule" 0. In this paper we present an alternative derivation of the gauge 
invariant path integral by applying a gauge covariant method developed long ago by 
Schwinger to study QED with background fields ||. In this method instead of the 
"x — p" ordering we use a time ordering procedure to obtain the configuration space 
path integral. 

The Hamiltonian operator for a charged particle in an arbitrary stationary magnetic 
field is 

H=—, (1) 
2m' V ; 

with Ti' L = p 1 — eA 1 (i=l,2,3) obeying the algebra (h = 1) 



[ki, TTj] = ie(diAj - djAi) 



(2) 



and the transition amplitude between the position eigenstates is given by 



(x", T|x', 0) = (x"|e 



iHT 




(3) 
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One way to write a path integral for (x",T|x',0), is to decompose it according to 
Dirac ||] as 



(x",T|x',0) = / dx N ...J dx 1 (x",T\x N ,t N )(x N ,t N \x N _ 1 ,t N _ 1 ) . . . (x 1 ,t 1 \x',0) . (4) 

Then we take T = (N + l)e and t n = ne (n = 1, . . . , N) letting N — > oo and s — > 
with T fixed 0. Next we obtain the small time propagator (x n+1 , t n + e|x n , £ n ), insert 
it in (|J) and try to perform the infinite dimensional integral. 

To evaluate the small time propagator we use a method developed by Schwinger in 
his early investigations on effective actions ||. We now write the transition amplitude 
as (x",T|x',0) = exp(iW), where W(x",x';T) is a complex function of the end point 
coordinates and time. Defining the expectation value of an observable O by 

(O)^(x",T|d|x',0)/(x",T|x',0), (5) 

it is easy to verify that W is determined by the following equations 



dW(x",x'; 


T) 


dT 




dW(x",x'; 


T) 


dx'( 




dW(x",x'; 


T) 


dx'i 




W(x",x' 


;0) 



= (H(x(T)MT))), (6) 
= (HT)) + eMx>>), (7) 

: -(7r i (0))-eA J (x / ), (8) 
= -iln5 3 (x" -x'). (9) 

To solve this problem Schwinger noticed that the above equations relate the tran- 
sition amplitude to the solution of the Heisenberg equations for x(T) and it(T). If we 
solve for tt{T) in terms of x(T) and x(0) and insert this, in a time ordered fashion, on 
(§)~(@) we are l e ft with a set of first order equations to integrate. For x(e) we have up 
to second order in e 

xHe) = e iflE xH<S)e- i6E ~ x*(0) + 7 ^-e + -^(2F lk n k (<d) - id k F ik )e 2 , (10) 

m Am 1 



with F tk = F* fc (x(0)). Inverting the above equation to get 7r 4 (0) in terms of x(e) and 
4(0) 

(x l {e) -x l (0)) e 
e 2 



11) 



using the fact that ((x(e) — x(0)) 2 ) is of order e ||, we see that if we take (7r(0)) the 
terms in the above expansion are respectively of order 1 / and yfe. The second term 
although small in comparison with the first, will give a relevant contribution when used 
in (|7|) and (§) to evaluate iy(x n+ i, x n ; e). From 7r(0) we have by time evolution: 



7T (e) ~ m 



- x 4 (0)) , e_,, . 4/n ,,- it 
v — + ~(ar(e) - x fe (0))F Jfe . 



(12) 
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Using the above expression for 7r(0) or 7r(e) in if, in a time ordered manner, we are 
ready to integrate (Q) 

^ 2 3 



W (x n+ i,x n ;e) ~ — h«-ine + <P(x n+ i,x n j , 



(13) 



where we used that [x l (e), Xi(0)] ~ —3ie/m and $ is a time independent function of 
the end point coordinates x n+i and x n . Inserting the above W in (|7|) and (H) we have 
for $ (remember (x n+ i — x n ) 2 ~ e) 



<9$ 



with the solution 



•''//-!- 1 - )' 7M f X; " :) + ^(x n+1 ) + ^(x r 



<9x 



71+1 

H x n) 



+ o( £ ), 

+ 0(e), 



(14) 
(15) 



$(x n+ i,x n ) = -(x n+ i - x n ) • (A(x n+ i) + A(x n )) + C. 



where C is a constant determined by 



C = -illn{— ). 
2 y 2m 



Finally, the small time propagator is 



m 
2ixie 



is 



(16) 



;i?) 



(18) 



3 



with 

S = m(Xn+ 2 X ~ Xn)2 + |(x n+1 - x n ) • (A(x n+1 ) + A(x n )) , (19) 

that is precisely the midpoint expansion rule. 

We see that by applying a gauge covariant method we were able to obtain the the 
"midpoint rule" in a natural way In our deduction two facts were essential, first the 
Brownian like relation (x n+ i — x n ) 2 ~ e and second the gauge covariance conditions 
(0) and (§). The authors are grateful to the CNPq (Brazilian Research Council) for 
the financial support. 
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